Dynamics of Conformal Maps for a Class of Non-Laplacian Growth Phenomena 
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Time-dependent conformal maps are used to model a class of growth phenomena limited by 
coupled non-Laplacian transport processes, such as nonlinear diffusion, advection, and electro- 
migration. Both continuous and stochastic dynamics are described by generalizing conformal- 
mapping techniques for viscous fingering and diffusion-limited aggregation, respectively. A gen- 
eral notion of time in stochastic growth is also introduced. The theory is applied to simulations of 
advection-diffusion-limited aggregation in a background potential flow. A universal crossover in mor- 
phology is observed from diffusion-limited to advection-limited fractal patterns with an associated 
crossover in the growth rate, controlled by a time-dependent effective Peclet number. Remarkably, 
the fractal dimension is not affected by advection, in spite of dramatic increases in anisotropy and 
growth rate, due to the persistence of diffusion limitation at small scales. 



Laplacian growth models describe some of the best 
known phenomena of pattern formation far from equi- 
librium, including continuous dynamics such as viscous 
fingering Q and (quasi-static) dendritic solidification 
and stochastic processes such as diffusion-limited aggre- 
gation (DLA) [j| and dielectric breakdown In this 
class of models, the interfacial velocity is determined 
by the normal derivative of a harmonic function, so the 
powerful technique of conformal mapping has been used 
extensively in two dimensions. Time-dependent confor- 
mal maps are used in the classical analysis of continuous 
Laplacian growth 0, 0, 0} ) and the analogous method of 
iterated conformal maps has recently been developed for 
stochastic Laplacian growth [a, y, llflj . 

In spite of the broad relevance of these models, real 
growth phenomena often involve non-Laplacian transport 
processes, such as advection or electro-migration coupled 
to diffusion 0, . Much less theoretical work exists in 
such cases, in part because conformal mapping would ap- 
pear to be of little use for non-harmonic functions. An 
exception is the recent use of streamline coordinates for 
dendritic solidification in a potential flow, but it turns 
out that advection has no effect on the shape of an in- 
finite dendrite ^jjj]. Stochastic conformal-map dynam- 
ics, however, has not yet been formulated for any non- 
Laplacian transport process (although iterated conformal 
maps have been used in a recent model of brittle fracture 
with a bi-harmonic elastic potential (l4jl- 

In this Letter, we formulate the dynamics of conformal 
maps for growth phenomena limited by non-Laplacian 
transport processes in a recently identified conformally 
invariant class [l5j|. We consider both continuous and 
stochastic growth from a finite seed to allow nontrivial 
competition between different transport processes. To 
illustrate the theory, we study fractal growth driven by 
advection-diffusion in a potential flow. 

Transport Processes — Consider a set of scalar 'fields', 
4> = {4>i, <f>2, ■ ■ ■ , 4>M}t whose gradients produce quasi- 



static, conserved 'flux densities', 

M 

F i = 5^C , «V0 J) V-F, = (1) 

i=i 

where the coefficients, {CV,- (</>)}, may be nonlinear func- 
tions of the fields. This general system contains a number 
of physical cases [l5j: (M = 1) simple nonlinear diffusion, 
Fi = —D(c) Vc, where c is a temperature or particle con- 
centration and -D(c) the diffusivity; (M — 2) advection- 
diffusion, Fi = -D(c)Vc + cV</>, of a scalar c in a po- 
tential flow, u = F2 = V</>, as well as ionic transport in 
a supporting electrolyte; and (Af > 2) ionic transport in 
a neutral electrolyte, Fj = — £)j(cj)Vcj — bi{ci)qiCi^ '<f>, 
where c$, Di, hi, and qi are respectively the concentra- 
tion, diffusivity, mobility, and charge of ion i, and <fr is the 
(non-harmonic) electrostatic potential, determined im- 
plicitly by electro-neutrality, J2i=i 1i c i = 0- 

In planar geometries, it is convenient to represent a 
vector, F = (F x ,F y ), as a complex scalar, F — F x + iF y , 
so Equation (JTJ takes the form, 

M 

in the z = x + iy plane, where V = + i-j^. Under a 
conformal mapping, w = f(z), Fi transforms as 

F i (z,z)=J 7 (z)F i {w,w). (3) 

(like V 0), and Re{VFi) = is unchanged. Therefore, 
even though the solutions (depending on z and z — x—iy) 
are not harmonic functions, the usual trick of conformal 
mapping still works 15]: If <p(w, w) solves Eq. (J2J in one 
domain, Sl w , then <f>(f(z),f(z)) solves Eq. in another 
domain, fl z = f^ 1 (fl w ) (with appropriately transformed 
boundary conditions). 

Interfacial dynamics in the plane can be elegantly de- 
scribed by a conformal map, z = g(w,t), from Q w , the 
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exterior of the unit circle, to f2 z (i),to the exterior of the 
(singly connected) growing object [a Is I& Ui| • Since 
the map must be univalent (1-1), it has a Laurent series, 

g(w,t) = A 1 (t)w + A (t) + ^^- + ..., \w\>l, (4) 

w 

where Ay(t) is real and defines an effective diameter of 
&z(t) @ • As described above, the fields satisfying (JTJ in 
Cl z (t) are easily obtained from the inverse, w = f(z,t), 
once the same equations are solved in Q w . As in Lapla- 
cian growth, the removal of geometrical complexity from 
the transport problem is a tremendous simplification. 

Boundary Conditions (BC) Along the moving 
boundary, <9f2 z (i) we consider generalizations of Dirichlct 
(<j>i = 0) and Neumann (n • = 0) BC [3: 



Ri (<f>(z,zj) = or h-Fi = Re [n{z)F i {z,z) y j = 0, (5) 

respectively, for z € dfl z (t), where n = n x + wi y rep- 
resents the outward normal, n, and Ri{4>) is a func- 
tion of the fields. The former BC express interfacial 
equilibrium for 'fast reactions' (compared to transport 
rates) , while the latter expresses impermeability to flows 
(F = u = V(/>) or flux densities of non-reacting species. 
Due to Eq. Q, these BC are the same, 

Ri(<f>{w,w)) =0 or n-Fi=Re(wFi{w,w)) = 0, (6) 

for w £ dd w , since n(w) — w — n(z)f'(z)/\f'(z)\. 

Far from the growth, we assume either constant values 
of the fields (e.g. temperature, concentration) or given 
flux (or flow) profiles which drive the growth: 

or Fi{z,z)r.Fr{z,z) (7) 

as \z\ — > oo (where F?° could also vary in time). The for- 
mer BC also remain the same after conformal mapping, 
but the latter is transformed by Eqs. © and Q: 

&(tu,tD)->#° or Fi{w,w) ~ TxF^iAxW^A^) (8) 

as | to | — > oo. Through Ay{t) (= Ay(t)), the fields and 
fluxes in fl w vary with the diameter of the growth, Q z (t) . 

Continuous Dynamics — Suppose that a Lagrangian 
boundary point, z{t) G dVL z (t), moves in the normal di- 
rection with (complex) velocity, 

M 

v = z t = ancr, a — Re(nQ)), Q = (9) 

i=l 

where Q is a flux density causing growth, a is a constant, 
and Bi(4>) may be functions of the fields. This generalizes 
Stefan's law (v = aft- V</>), e.g. to electrodeposition |T^ 
(where Q is for the depositing ion). 

The classical analysis of viscous fingering is easily gen- 
eralized to Eqs. @, ©, 0, and @. Let wit) be a 
'marker' for z(t) — g(w(t),t) on ft w 7]. Substituting 



Zt = g' w t + gt into Eq. ©, multiplying by wg 1 , taking 
the real part, and using Ke(ww t ) = for |u>| 2 = ww = 1, 
we arrive at an evolution equation for the conformal map, 

Re(w7g 7 5t ) = a a(w,Ai(t)) for \w\ = 1 (10) 

where o~(w, Ay (i)) = Ke(wQ(w, W)) is the normal flux 
density on Q w , which depends on Ay(t) through Eq. JHJ. 

The evolution equation for radial Laplacian growth 
(e.g. viscous fingering) 0,0, Ke(wg'gt) = 1, corre- 
sponds to the special case of uniform flux, a = constant. 
This dynamics is known to preserve number of pole-like 
singularities (inside the unit circle) for a wide class of ini- 
tial maps, g(w, 0). Except for some elementary maps, e.g. 
circles and ellipses, which preserve their shapes, arbitrary 
smooth initial interfaces develop singularities (cusps) in 
finite time |(J. 

In our generalized models, the time-dependent, non- 
uniform flux, er, in Eq. l|lUfl changes the analytic struc- 
ture of an initial map (e.g. number of poles), so even 
circles and ellipses become distorted. This raises interest- 
ing questions about finite-time singularities, e.g. What is 
the fate of solidification from a circular seed in a flowing 
melt (with a described below)? Does advection gener- 
ally enhance or retard the formation of singularities? We 
leave these questions for future work and focus here on 
non-Laplacian fractal growth. 

Stochastic Dynamics — Suppose that the domain, 
f^z(ijv), grows from its initial shape, Cl z (0), at times 
ti,t%, ■ ■ . , ijy by discrete 'bumps' representing N parti- 
cles of characteristic area, A |a @. Since our mod- 
els exhibit non-trivial time-dependence (see below), we 
introduce time into the usual morphological model by 
replacing Eq. Q with p{z, fjv) = aa(z) /A , for z S 
dn z (tN-i), where p{z,tj^) \dz\dt is the probability that 
the iVth growth event occurs in the boundary element 
(z,z + dz) in the time interval (t,t + dt). The waiting 
time, ijv — ijv— l; is then an exponential random variable 
with inverse mean, 

^ = T- / *( z ) \ dz \ = T- 1^ °{z l6 :MtN-y)) d9, 
*o Jan,(t N _i) Jo 

(11) 

where we use \dz\ = \dw\/\f'\ and a(z) = \f'(z)\a(w) 
from Eq. J2J to transform to d£l w where w = e and 
\dw\ = d9. The probability that the position of the Ath 
growth occurs in (z, z + dz) C d£l z {tM-y) is, pN(z)\dz\ = 
f- o r N a{z)\dz\ = P N {d)d6, where 

PN(6) = ^T N a{e ie ,Ay(t N - 1 )), e l6 e Q w (12) 

is the probability measure for angles on the unit circle. In 
DLA [8j, pn(z) is the harmonic measure; Pn{@) = (27r) _1 
is uniform; and tn is not defined. 

It is now straightforward to generalize the Hastings- 
Levitov DLA algorithm to our non-Laplacian models. 
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FIG. 1: Evolution of the flow (yellow streamlines) and concentration (color contour plot) in £l m (top) and in Q z {t) (bottom) 
during advection-diffusion-limited aggregation of N = 10,000 particles (A = l? = 1) with Pe D = 0.05. The columns from left 
to right correspond to Pe = 0.1, Pe = 1, and Pe = 10. The lower plots are shown rescaled by A\{t$) = 2, A\{t\§z) = 20, and 
Ai(t962i) = 200, the effective diameters of the fractal cluster. 



The univalent map z = g{w, ijv) from Q w to f2 z (ijv), 
the exterior of an TV-particle aggregate, is constructed 
iteratively from a two-parameter basic map that places a 
'bump' of area Ajv around an angle On on the unit circle. 
In order to fix the area A D of a new bump on dil z (tN) 7 
the size of its pre- image, An, on dtt w is divided by the 
Jacobian of the previous map, Afc = \ / \g' (e l9 N , ijv— 1)| 2 - 
The first difference with DLA is that the angle, 9jsr, is 
chosen according to the time-dependent (non-harmonic) 
measure, Pjy(6), in Eq. I|12l) . The second difference is the 
evolving waiting time, tn, in Eq. (|llfl . 

Advection- Diffusion- Limited Aggregation — (ADLA) 
As an example, we consider the stochastic aggregation of 
particles around a circular seed of radius, L Q , limited by 
advection-diffusion in a uniform potential flow of speed 
U and concentration C. In £l z (t) the transport problem 
has the usual dimensionless form, 



Pe V(/> • Vc = V 2 c, V 2 </> = 0, 
c = 0, n ■ V</> = 0, a = n ■ Vc, 
c — * 1, V</> — * x, 



z€Sl 9 (t) (13) 
z e dfl z (t) (14) 
|z|^oo (15) 



where x, 0, c, and a are in units of L a , UL D , C, and 
DC/L , respectively, and Pe D = UL a /D is the initial 
Peclet number. In complex notation, we solve Eq. @ in 
n z (t) with BC © and Q for the 'fluxes', F x = -Pe Vc+ 
cV(j) and F2 — u = V0. In fl W) Equations J5J) and © 
are the same, but the BC calls for a background flow 
speed at |tu| — * 00 that diverges with A\(t). 

It is natural then to rescale the w-velocity by A\ (t) to 
fix the background speed at unity, and instead solve the 



same equations in £l w as in £l z (t) with a time-dependent 
Peclet number, Pe(t) — Ve A\(t). Since A\(t) is an ef- 
fective diameter for Q z (t), the theory has shown us how 
to properly define Pe (which is not obvious for fractals). 
We also see that advection eventually dominates diffusion 
since Pe(i) — * 00, so we expect to see a crossover from 
DLA to new advection-limited dynamics, as shown in the 
simulation of Fig. ^ which we now discuss in detail. 

With the scalings above, the velocity potential in fl w 
has the usual harmonic form, 4>[w) — Re(ui + 1/w), for 
potential flow past a cylinder, but the non- harmonic con- 
centration, c(w,w, Pe(t)), cannot be expressed in terms 
of elementary functions. However, asymptotic expan- 
sions (for fixed w) can be derived \u\ . 



c(w, w, Pe) 



-Re (log w/ logPe) 



rf 



v^Im 



Pe < 1 
Pe > 1 



(16) 

The low-Pe approximation, the familiar harmonic field 
of Laplacian growth, is valid out to a 'boundary layer at 
00', while the high-Pe approximation is valid away from 
a thin wake on the positive real axis (the branch cut for 
[3. From Eqs. (JTTJ, and JTHJ we also obtain, 



tn 



— log Pe 



Pn(0) 



l/2n 
1 c ,;„ ' 



Pe < 1 
Pe > 1 



(17) 



where t/v is measured in units of (X /D)/(aC). Even 
for Pe G -C 1, the ADLA dynamics smoothly crosses over 
from the DLA 'unstable fixed point', Pe(t) <C 1, to a 
new advection-dominated 'stable fixed point', Pe(t) ^> 1, 
described by Eq. JT7J . 
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FIG. 2: Log-log plot of the Laurent coefficients, A\(t) 
and \Ao (t) | versus N(t), averaged over 30 ADLA simulations 
(Pe D = 0.05, A = L„ = 1), compared with analogous results 
for DLA without advection. 

To investigate this crossover, we accurately obtain the 
normal flux, o~(w, Pe (£)), on d£l w by interpolating static 
numerical solutions for a range Pe. Following the meth- 
ods above, we then perform many ADLA simulations of 
N > 10 4 particles for various Pe Q and \ jl? . Details will 
be given elsewhere, but here we briefly discuss morpho- 
logical changes and growth rate. 

The Laurent coefficients in Eq. contain morpholog- 
ical information |3, |9( ■ As illustrated in Fig. [21 for all 
Pe , the diameter of the cluster, Ai(ijv) ~ N x / D > ', re- 
markably maintains the same fractal dimension as DLA, 
Df = 1.71, for all Pe(t) and N(t) in the scaling regime, 
Ai(tN) S> 1. (As a check, we obtain the same Df from 
the radius of gyration.) A posteriori, this can be under- 
stood by noting that the stable fixed point has a growth 
measure, Pn(@) oc sin6*/2, which is differentiable, and 
thus locally constant (as in DLA), everywhere except at 
the rear stagnation point, 9 — 0. Physically, this simply 
means that diffusion always dominates advection at small 
scales. We conjecture that the same Df holds for any 
non-Laplacian dynamics in our class, if ImiN^oo Pn(9) 
is continuous and almost everywhere differentiable. The 
surprising universality of Df = 1.71 makes its exact value 
seem quite fundamental. 

Of course, as seen in Fig. ^ the growth is highly 
anisotropic and moves toward the flow for Pe(t) ^> 1. 
This is seen in scaling of the next Laurent coefficient, 
A (t), the 'center of charge' 0, 

( N a Pe «: 1 

\Ao(t N )\~[ Nl/Df (18) 

which crosses over from DLA scaling (2a = 0.7 |3|) to the 
same scaling as Aiit^), as shown in Fig.El Elsewhere, we 
will show that |Ao(t)|/Ai(t) tends to a universal function 
of Pe(i) after initial transients vanish, Ai(t) 1, e.g. 
lim t ^ oo |4>(t)IMi(t) = 0.6. 

The expected total mass versus time, N(t), also un- 
dergoes a crossover. Using Eqs. ifTTJl and (fT7| and in- 
tegrating dN/dt = r w 1 (Pe(t)) = r-^PeoiV 1 /- /) yields, 
N(t) oc t for t < 1 and N(t) oc ^d / /(2D / -i) for f > ^ 



where again only Df matters. The cluster diameter Ai(t) 
switches from t 1 / 01 to t 2 /( 2D /~ 1 ' scaling in time, even 
though the scaling with N(t) does not change. One 
should also bear in mind that tn — > as Pe(i) — > oo, 
so the quasi-static, discrete-growth approximation must 
eventually break down (although this is delayed in the 
dilute limit, aC « 1). 

In summary, we have formulated conformal-map dy- 
namics for a class of planar growth phenomena limited 
by non-Laplacian transport processes. Although various 
simplifying assumptions were made, the method enables 
very efficient simulations of fractal-growth phenomena, 
such as ADLA, which it seems could not be achieved by 
more efficient way. By describing competing transport 
processes, it also paves the way for new studies of gen- 
eral crossover phenomena in pattern formation. 
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